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Wind  velocity  parameters  in  the  planetary  boundary  layer  are  measured  using 
spatial  inhomooeneities  in  the  natural  aerosol  content  as  tracers.  Lidar 
measurements  of  these  spatial  inhomogeneities  are  cornlated  in  time  and 
space  to  remotely  estimate  the  mean  speed,  direction  and  the  rms  speed  variant 
of  the  wind  velocity.  Aerosol  inhomogeneities  with  horizontal  dimensions  of 
60m-0.1km  were  found  to  be  efficient  tracers  of  the  wind.  The  coherence  be¬ 
tween  spatially  and  temporally  separated  lidar  profiles  of  these  spatial _ ^ 


aerosol  Inhomooeneltles  Is  used  to  determine  the  perpendicular  components  alono 
and  across  the  beam  pattern  of  the  horizontal  wind  velocity--  The  decrease  in 
coherence  with  increasino  time  and  spatial  separation  is  used  to  estimate  the 
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20.  ABSTRACT  (cont) 


rms  windspeed  over  the  wavenwnber  reoion  .001  '  -  .07(m  ).  Comparisons  be¬ 
tween  lidar  horizontal  wind  measurements  with  tower  ai. vnomcters  and  pilot  bal¬ 
loons  show  that  lidar  can  remotely  sense  the  boundary  layer  flow. 
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1.0  Introduction 

Optical  remote  sensing  of  the  boundary  layer  flow  can 
potentially  provide  the  spatial  and  temporal  detail  needed 
in  meteorological  application  in  this  region.  Convective 
boundary  layer  studies,  air  pollution  monitoring,  and  avia¬ 
tion  operations  require  knowledge  of  the  wind  velocity 
variability  and  atmospheric  dispersal  mechanisms  to  gain  a 
better  understanding  of  the  physical  principles  acting  in 
the  boundary  layer.  The  application  of  laser  technology  to 
measure  the  wind  velocity  can  provide  wind  information  with 
greater  sampling  statistics  and  negligible  influence  on  the 
flow  pattern  over  conventional  wind  measurements. 

A  number  of  different  techniques  have  been  investigated 
for  remote  determination  of  wind  velocities  using  optical 
methods.  Heterodyne  determination  of  the  Doppler  shift  in 
the  received  scattered  light  from  transmitted  laser  pulses 
have  been  used  by  Lawrence  et  al.  (1972)  and  Benedetti- 
Michelangeli  et  al.  (1972)  to  calculate  the  wind  velocity 
component  along  the  laser  path.  Kjelass  and  Ochs  (1974) 
determined  the  mean  horizontal  velocity  and  divergence  near 
the  surface  over  a  300  m  equilateral  triangle  by  correla¬ 
tion  of  scintillation  patterns.  Xuev  et  al.  (1977)  have 
used  the  correlation  of  two  spatially  separated  laser  beasi 
signals  to  determine  scale  sizes  and  life-times  of  aerosol 


Inhonoceneities.  Correlation  of  lidar  return  si^nnls 
scattered  from  aerosol  density  inhomogeneities  have  been 
used  to  estimate  mean  horizontal  wind  velocities.  Armstrong 
et  al.  (1976)  have  obtained  average  wind  measurements  over 
a  time  scale  of  seconds  from  a  distance  of  2^0  m. 

Eloranta  et  al.  (1975)  have  used  this  procedure  to  estimate 
the  average  wind  velocity  at  distances  of  several  kilo¬ 
meters.  Lidar  returns  for  large  scale  (60  m  -  1.0  kn) 
spatial  Inhomogeneities  in  naturally  occurring  aerosol 
content  are  used  in  this  study  to  determine  two  horizontal 
wind  velocity  components  and  also  an  estimate  of  the  rms 
wind  speed. 

Fig.  1  represents  a  schematic  diagram  of  the  procediure 
used  in  obtaining  lidar  measurements  of  the  spatial  aero¬ 
sol  inhonogeneities  at  different  times  and  azimuth  direc¬ 
tions.  Ihcsc  spatial  aerosol  i nhomogenei t ics  are  primarily 
due  to  nonunifom  aerosol  sources  and  turbulent  mixing  of 
particulate  matter  into  the  atmosphere.  Motion  of  the 
wind  can  be  inferred  by  detecting  the  movements  of  these 
aerosol  inhomogt neities  as  they  are  advected  through  one 
lidar  sample  volume  to  the  next.  Previous  studiea  by 
Eloranta  et  al.  (1975)  end  Leuthner  (1976)  have  shown  that 
velocity  information  can  be  obtained  by  correlation  between 
two  lidar  radial  aeroaol  inhomogeneity  profiles  separated 
in  time  and  space.  They  calculated  a  two-dimensional 


i 


?igur«  It  Sohcmatie  diagram  of  lidar  wind  maasuramenta 
using  natural  aarosol  as  traears*  Diagrams  ara  at  tiiM 
t  and  t-f-dt*  lidar  baam  has  morad  in  azimuth  diraetion  in 
this  tisia.  Aarosol  motion  is  astimatad  by  eorralating 
tha  lidar  inhomoganaity  profilas  batuaan  two  diffarant 
baams»  ^  $  t  and  ^  ,  t-t-at. 
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lag'Croaa  correlation  coefricient  in  time  lag  ^t  and  radial 
apatial  lag  from  lidar  returna  at  different  azimuth 
anglea  aa  ahown  in  Pig.  1.  The  maximum  in  the  lag^croaa 
correlation  coefficient  was  used  to  estimate  the  mean  hori« 
zontal  wind  velocity.  Because  turbulent  nixing  changes  the 
spatial  distribution  of  individual  aerosol  inhonogeneitiea* 
velocities  calculated  in  this  manner  will  be  biased  toward 
higher  values*  see  Briggs  et  al.  (19^0).  Kunkel  (1976) 
modeled  the  spatial  aerosol  inhomogeneities  as  Gaussian 
structures  being  advected  by  a  Gaussian  velocity  distribution. 
This  model  estimates  the  mean  velocity  and  the  wind  speed 
variance  from  the  lag-cross  correlation  coefficient  values  for 
different  lateral  separations  using  the  method  developed  by 
Leuthner  (197b). 

The  coherence  which  describes  the  correlation  at 
various  radial  wavenumbers  k  is  used  in  this  study  for 
velocity  calculations.  Spectral  radial  velocity  components 
are  calculated  from  the  spatial  phase  difference  at  each 
spectral  wavenumber  between  inhomogeneity  profile a  due  to 
their  radial  drift.  The  lateral  spectral  velocity  components 
are  calculated  from  the  time  lag  for  each  wavenuaiber 
component  to  drift  fr<»  one  azimuth  directicn  to  the  next* 
of.  Pig.  1.  The  mean  value  of  each  velocity  component  is 
obtained  by  averaging  over  wavenumbera.  k  model  aimilar  to 
that  of  Kunkol  (1976)  is  used  to  correct  for  the  effects  of 


turbulence  end  estimate  the  rms  wind  speed. 

An  increase  in  computational  efriciency  is  obtained  in 
transforming  to  radial  wavenumbers  k  by  using  Fast  Fourier 
Transform  (FKT)  in  calculation  of  spectra  and  coherence  of 
spatial  aerosol  inhomogeneity  profiles.  The  FFT  computes 
m  discrete  wavenximbers  which  depend  on  the  data  segment 
length  and  point  separation.  (For  this  study  k^  «  ^/960  m* 
i  3  1,  32.)  This  method  can  be  5  to  30  times  more  efficient 
computationally  over  the  Blackman-Tukey  method,  see  Cooley  et 
al.  (1970).  Filtering  the  spatial  series  to  reduce  the 
effects  of  noise  is  also  more  efficient  in  wavenumber  space. 

A  method  has  been  developed  to  optimise  the  filtering  of  the 
spatial  series  to  maximise  the  signal  to  noise  ratio  in 
changing  atmospheric  conditions. 

This  study  has  been  organised  in  the  following  manner. 
The  theoretical  procedure  used  to  determine  the  wind 
velocity  parameters  is  presented  in  Chapter  2.  Derivations 
and  computational  details  are  left  to  the  appendices  for 
thoae  with  further  interest.  Chapter  3  describes  the 
instrumentation  used  during  this  study.  The  re  axil  t  a  of 
experiisenta  conducted  at  Madison,  Wisconsin  and  the  Whits 
Sands  Missile  Range,  New  Mexico  are  shown  in  Chapter  U. 

The  conclusions  are  stated  in  Chapter  5* 
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2.0  Theor«tieal  Method  for  Wind  Velocity  Measurement 

This  chapter  presents  the  theoretical  method  used  in 
this  thesis  for  remote  horizontal  wind  velocity  measurements 
using  a  monostatie  lidar.  It  is  divided  into  four  sections} 
the  procedure  for  obtaining  measurements  and  processes 
required  to  calculate  coherence  values  and  spectral  velocities, 
estimation  of  the  spectral  signal  to  noise  ratio  and  optimal 
filtering, wind  velocity  parameters  calculated  from  the 
correction  for  turbulence,  and  uncertaluties  in  lidar  wind 
measurements.  Mathematical  details  and  derivations  are 
contained  in  the  appendices. 

2.1  Procedure}  Obtainittg  Lidar  Wind  Measurements 

The  procedure  to  obtain  measurements  used  to  infer 
horizontal  wind  velocity  from  the  notion  of  spatial  aerosol 
inhomogeneities  is  described  in  this  section.  A  method  is 
presented  to  calculate  these  spatial  inhomogeneities  in  the 
mean  aerosol  content  from  lidar  measurements.  The  experi* 
mental  arrangement  for  obtaining  these  measurements  is  also 
shown.  Borisontal  wind  velocity  information  is  calculated 
from  tbs  eoherenoe  between  lidar  measured  radial  profiles  of 
spatial  aerosol  inhomogeneities. 

Lidar  is  an  active  remote  sensing  device  used  to  probe 
the  atmosphere.  A  laser  pulse  is  transmitted  into  the 
atmosphere  and  measurements  are  made  of  the  photons  scattered 


through  lUO  decrees  by  atmospheric  constituents  into  a 
receiver  whose  axis  is  aligned  parallel  to  the  laser  axis. 
These  detected  scattered  photons  are  converted  into  an 
electrical  signal  which  is  proportional  to  the  amount  of 
scattered  laser  light.  The  equation  for  singly  scattered 
monostatie  lidar  return  power  is^ 


€  C 


fpct  jyj 
y  7r 


^  .  • 


2.1.1 


where  R  is  the  radial  distance  along  the  propagation  path 
R*  is  the  integration  variable  along  R 
P^(R)  is  the  instantaneous  received  power  from  range  R 
ia  the  tranamitted  laser  pulse  energy 
Ap  is  the  area  of  the  receiver  telescope 
e  is  the  speed  of  light 

>  is  the  scattering  oross*aeetion  per  unit  volume 
/3g<^^  ia  the  extinction  oross*aection  per  unit  volume 
the  backscatter  phase  tunction 
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A  modified  form  of  the  lidar  equation  presented  by 
Coilia  (19691. 
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A  description  of  the  University  of  Wisconsin  lidar 
system  is  presented  in  section  3«1*  Lidnr  returns  are  pre- 
processed  by  this  system  in  the  following  -manner.  Loga¬ 
rithmic  amplification  compresses  the  dynamic  range  of  the 
signal.  The  signal  is  digitized  into  10  bit  words  at  a 
lOLHz  rate  which  yields  a  15  m  resolution  along  the  laser 
beam  path,  biuar  returns  are  correcuea  lui-  tne  “range 
square  attenuation”  and  normalized  by  the  transmitted  laser 
pulse  energy  using  a  PDF  11/UO  minicomputer.  Lidar  returns 
analyzed  in  this  study  are  given  by  K  the 

natural  logarithm  of  the  quantity  (Pp(K)  R^/Eq); 


2.1.2 


where  t  is  the  time  at  which  the  nth  laser  shot  occurred 

D 

is  the  range  to  the  1^^  data  point 
(t^»  R^)  is  the  lidar  return  analyzed  at  azimuth 
angle  / 

K  is  a  systems  parameter  constant 

A.J  volume  backscatter  coefficient 

1  TT 


Spatial  inhomoceneitioa  of  the  atmospheric  aerosol 
content  can  be  obtained  from  lidar  measurements.  Spatial 
fluctuations  in  concentration,  sire  and  shape  of  aerosol 
particles  produce  spatial  variations  in  the  backscatter 
coefficient  .  The  attentuation  term  in  Eq.  2,1,2 

^  is  a  monotonically  increasing  function  with 
range  and  is  assumed  to  be  only  slowly  varying  in  time. 
Estimation  of  the  spatial  inhomogeneities  in  the  mean  aerosol 
distribution  can  be  obtained  by  calculating  radial  deviations 
from  a  time  centered  running  mean  radial  profile  of  Eq,  2,1,2, 
for  each  azimuth  angle.  Spatial  deviations  from  the  mean 
aerosol  profile  calculated  in  this  manner  remove  the  effects 
of  slow  temporal  variations  in  the  extinction  (cf,  Eq,  2,1.2,) 
in  correlation  between  profiles.  These  deviations  from  the 
mean  aerosol  profile  are  given  by 


where  is  the  number  of  lidar  profiles  used  to  calculate 
the  average  lidar  profile 


is  the  deviation  in  Yv  from  the  average 


profile 
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Tho  avcrapo  vnlue  oT  Y'^  M^)  alonp.  the  rah^c  is 

renoved  before  harnonic  analysis  to  reduce  tho  leakage  of 
zero  wavenunber  conponent  into  other  wavenurabor  bands. 
Radial  profiles  of  the  spatial  aerosol  inhonogeneities 
analyzed  in  this  thesis  are  defined  as 


»7  J 

:  e. 


N 


2.1.4 


Where  N  is  tho  nunber  of  data  points  along  the  radial  range 
If  is  the  first  data  point  of  this  segment. 

This  procedure  also  suppresses  spurious  harmonics  generated 
by  discontinuities  in  the  aerosol  structure. 

Lidar  measurements  of  the  spatial  distribution  in 
aerosol  content  are  obtained  from  a  sequential  three  angle 
azimuth  PPI  scan  at  a  fixed  elevation  angle,  see  Pig.  2. 

The  sequential  scanning  procedure  obtains  lidar  profiles  of 
the  aerosol  content  at  azimuth  angles  etc.. 

Short  data  segments  along  two  beam  paths  can  be  treated  as  if 
they  are  parallel  for  small  angular  separation,  , 

and  tho  change  in  the  lateral  separation  at  the  endpoints  is 
small  compared  to  the  average  lateral  separation.  The  lateral 
separation  distance  Ax  is  defined  as 


AX  -  2  R  sin  (®</2) 


2.1.i> 


Pigur*  2t  Vertical  and  horizontal  riewa  of  lidar  wind 
measurement  scheme.  Lidar  n^soreaents  are  obtained  at 
three  azimuth  angles  (^2*  ^2*^^  with  an  elsTation  angle. 

The  wind  velocity  has  components  perpendicular  and 
parallel  to  the  center  azimuth  (^2^  ^ 

horizontal  range  is  given  as  Rw*Reos(9). 


I- nr 
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The  height  Z  of  a  data  point  at  a  range  R  ia 


Z  =  R  sin  (9)  2.1.6 

where f»is  the  elevation  angle  of  the  lidar  system.  For  the 
geometry  shown  in  Fig.  2,  wind  measurements  at  different 
height  levels  correspond  to  correlations  of  lidar  returns 
at  various  range  values. 

Spatial  aerosol  inhomogeneities  in  the  average  lidar 
radial  profiles  calculated  from  Eqs.  2.1.3  >  2.1.U  are 
Fourier  analysed  to  calculated  power  spectral  and  cross 
spectral  estimates  of  the  spatial  aerosol  inhomogeneities. 

The  details  of  the  analysis  in  computing  the  smoothed  power 
spectral  and  cross  spectral  estimates  of  Y'y»*  described  in 
Appendix  A.  These  smoothed  spectral  estimates  are  used  in 
the  correlation  of  spatially  and  temporally  separated  lidar 
returns  to  deduce  the  spatial  aerosol  inhomogeneity  motion. 

The  coherence  is  a  measure  of  the  degree  of  correlation 
between  two  spatial  series  as  a  function  of  radial  wavenumber 
k.  Coherence  values  of  data  segments  centered  at  a  range  R 
are  calculated  for  various  time  lags»  dt  and  for  three 
lateral  separations,  ax  «  Om,  ax  «  2  R  sin((  )/2}  and 
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CoL  * 


•V,'*  ’ 


2.1.7 


where  »  denotes  conplex  conju>;ation 
k  la  the  radial  wavenunber 

at  Is  the  time  lag  between  radial  profiles  at  ^ 
d  X  is  the  average  lateral  separation  between  radial 
profiles  at 

is  the  smoothed  spectral  estimate  of  the 

radial  profiles  of  the  spatial  aerosol 

inhonogeneities, 

^>0  ^ 

S  )  is  the  smoothed  spectral  estimate  of  the 

radial  profiles  of  the  spatial  aerosol  inhomogeneities 

i>  the  smoothed  cross  spectral  estimate  between 
radial  profiles  of  the  spatial  aerosol  inhomogeneities 
measured  by  and  Y*^  .  Y*Jf^  and  Y*^  are 

separated  by  an  average  lateral  distance  ox  and 
time  lag  at 

me  process  of  smoothing  the  spectral  estimates  is 


important  in  ealctUating  reliable  measurements  of  the 
coherence.  Smoothing  the  spectral  estimates  decreases  the 
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probability  that  ranUo.’n  data  will  yield  a  ni^cn  coherence 
value.  The  probability  of  raiidon  data  having  a  given 
coherence  value  is  (cf.  Eq.  A, 2. 2) 

P  =  2.1.b 

where  P  is  the  probability 

p  is  the  coherence  value 

edf  is  the  equivalent  number  of  degrees  of  freedom 
The  probability  of  random  data  having  a  coherence  value  of 
1.0  is  100%  in  the  limiting  value  for  an  equivalent  number 
of  degrees  of  freedom  equal  to  2.  For  data  presented  in  this 
study  coherence  values  greater  than  .17  have  a  probability 
of  being  random. 

Estimation  of  the  horizontal  wind  velocity  can  be 
obtained  by  reliable  measurements  of  the  coherence  between 
lidar  radial  profiles  of  the  aerosol  inhomogeneities. 
Individual  spectral  components  are  assumed  to  be  propagating 
with  a  spectral  velocity  V(k).  The  relative  phase  difference 
in  the  smoothed  aerosol  inhomogeneity  cross  spectrum  is 
proportional  to  the  distance  over  wiiich  the  spectral  component 
has  drifted  radially. 

At  A  2.1.9 

•»7r 

where  Vj,(k)  is  the  spectral  velocity  component  along  the 
beam  path,  cf.  Pig.  2.  The  solution  for  the  radial  spectral 


velocity  is 


IS 


4  V/  <i> 

"  WjTlt  2.1.10 

The  lateral  spectral  velocity  conponent  u(k)  can  be 
estimated  from  the  coherence  calculations  between  lidar 
aerosol  Inhooogeneity  profiles  at  different  azimuth  angles, 
ef.  Pig.  2.  A  maximum  coherence  value  in  time  lag  dt 
determines  approximately  the  time  for  a  spectral  component 
to  drift  the  average  lateral  separation  distance  ^  x.  An 
approximate  lateral  spectral  velocity  component  u*(k)  is 
given  by 

u»(k)  »  fix/dt  (k)  2.1.11 

re 

where  dt^(k)  is  the  time  lag  value  for  maximum  spectral 
coherence.  Random  changes  in  the  spatial  aerosol  inhomo¬ 
geneity  profiles  due  to  turbulent  mixing  cause  the  degree 
of  correlation  (coherence)  to  decrease  with  increasing  4X 
and  4t,  ef.  Briggs  et  al.  (1950).  Section  2.3  presents  an 
approximation  to  estimate  the  mean  lateral  velocity  from  the 
coherence  decay. 
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2.2  Signal  Spectrum  to  Noise  Spectrum  Estimation  and 
Filtering 

This  section  describes  the  estimation  of  an  optimal 
linear  filter  based  upon  the  spectral  characteristics  of  the 
lidar  radial  aerosols  profiles.  A  method  to  calculate  this 
optimal  filter  from  lidar  coherence  measurements  with  no 
lateral  separation  is  also  developed.  The  derivation  of 
this  procedure  is  detailed  in  Appendix  B.  This  procedure 
allows  the  spectral  characteristics  of  the  data  being 
processed  to  dictate  the  filtering  characteristics.  Examples 
of  the  effectiveness  of  this  procedure  are  included  at  the 
end  of  this  section. 

A  number  of  sources  contribute  to  noise  in  lidar  mea* 
surements.  Background  light  and  statistical  fluctuations  in 
detecting  scattered  photon  are  sources  of  noiae  in  lidar 
measurements.  Electronics  and  digitisation  noise  also 
contribute  to  the  total  noise  in  lidar  returns.  Rapidlj 
evolving  small  aerosol  density  structures  are  also  part  of 
the  backgroxind  noise.  Spatial  filtering  oi  ^idar  radial 
profiles  is  necessary  to  reuuce  tnese  noise  cottbriOubxuns 
from  the  signal  of  the  spatial  aerosol  inhomogeneities. 

A  filter  that  reproduces  the  signal  from  the  aerosol 
inhomogeneitiea  with  a  minimum  in  error  can  be  calculated 
from  knowledge  of  the  spectral  characteristics  of  both  the 
signal  S|(k)  and  noise  S^lk)  contributions,  cf.  Weinstein 
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and  Zubakov  (1962).  The  Fourier  transform  K(K)  of  the 
optimal  linear  filter  h(x)  is  given  by,  cf.  Eq.  B.7. 

S  (K) 

•Hk)  =  - ^  2.2.1 

■Sg(k)  ♦  S'j^(k) 

This  optimal  transfer  function  U(k)  is  derived  for  mea« 
surenents  where  the  signal  and  noise  contributions  are  not 
correlated.  The  noise  contributions  to  lidar  returns  are 
assumed  to  correlate  only  with  the  same  return  and  any 
correlation  between  different  lidar  profiles  is  primarily 
due  to  signal  contributions.  Power  spectral  and  cross 
spectral  estimates  of  the  aerosol  inhomogeneity  signals  are 
sufficiently  smoothed  to  decrease  the  probability  of  random 
statistical  correlations  (see  discussion  in  section  2.1  p.l3  ). 

An  estimate  of  the  ratio  of  the  signal  spectrum  to  the 
noise  spectrum  can  be  obtained  from  coherence  measurements 
between  spatial  series  in  the  limit  as  the  time  separation 
approaches  zero.  The  derivation  of  this  procedure  for 
estimating  the  signal  spectrum  to  noise  spectrum  ratio  from 
ooherenoe  measurements  is  detailed  in  Appendix  B.  Coherence 
values  are  calculated  from  experimental  lidar  measurements 
of  aerosol  inhomogeneities  for  the  lateral  separation 
Hx  s  0,  see  Pig.  2.  These  coherence  values  are  fxinctions  of 
wavenxunber  and  time  lag  only.  They  are  observed  from  Pig.  7 
to  be  approximately  Oausaian  in  time  lag  At.  Extrapolation 
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of  coherence  moasurencnte  at  each  wavenumber  k  to  =  0 
is  estinated  by  a  least  squares  regression  to  a  Gaussian  of 
the  form 

Coh(k,0,£it)  =  AlkJ  .  exp(-Btk)' fit^)  2.2.2 

where  A(k)  and  5(k)  are  the  regression  parameters  at  each 
wavenumber.  The  extrapolated  coherence  value  for  =  0  is 
just  the  amplitude  factor  A(x)  (note  that  A(k>  =  1  for 
noise  free  measurements).  Tne  ratio  of  the  signal  spectrum 
to  the  noise  spectrum  can  be  estimated  by 


S,(k) 


U  ‘‘  (»)  -  i)"' 


2.2.3 


where  this  equation  is  derived  in  appendix  B  as  Eq.  B.b. 
The  transfer  function  of  the  optimal  linear  filter  can  be 
calculated  from  knowledge  of  the  signal  spectrum  to  noise 
spectrxim  ratiOf  cf.  Eqs.  B.^  and  B.7. 


H(k)  =  (Coh(k,0,  =  A**  (k)  2.2U. 

The  filter  corresponding  to  this  transfer  function  will 
optimize  the  recovery  of  signal  information  of  spatial 
aerosol  inhomogeneities  from  lidar  aerosol  profiles. 

An  estimate  of  the  ratio  of  the  total  filtered  signal 
variance  to  noise  variance  (Sh'B)  is  calculated  here  by 
replacing  the  integrals  in  Eq.  B.IO  by  sums  over  all 


wavenur.bera 

h 

Sf.H  =  ^  j  2.2.i> 

Z  H  f&) 

4 

Where  the  sums  are  over  the  K  discrete  wavenamoers  obtained 
from  the  Fast  Fourier  Transform  (FPf). 

Figs.  snow  experimental  results  from  lidar  data 

obtained  on  January  19»  19Yb  at  White  Sands  Missile  Range, 

Kew  Mexico.  Throe  cases  are  shown  for  different  values  of 
SNR  to  illustrate  the  effects  of  noise  upon  the  data,  'llie 
cases  of  high  iSf«R=2.90),  intermediate  (SN’R=1.19)  and  low 
(SNR=:0.b5J  signal  to  noise  ratio  nearly  encompass  the  total 
range  observed  in  the  data  presented  in  this  thesis. 

The  power  spectra  of  the  spatial  aerosol  inhonogeneities 
calculated  from  Eq.  2.1.4  are  shown  in  Fig.  3.  These  spectra 
are  not  expected  to  follow  the  >5/3  power  law  since  the 
logarithm  of  the  lidar  returns  are  used  in  analysis  for 
increased  computational  efficiency  and  to  remove  the  effect 
of  slowly  varying  extinction. 

Pigs.  4  OAd  3  ere  similar  graphs  for  the  estimated 
signal  spectrum  to  noise  spectr\in  ratio  (Eq.  2.2.3)  end 
the  optimal  transfer  function  (Eq.  2.2.4)  for  the  data 
presented  in  Fig.  3.  Pig.  U  shows  that  larger  values  of 
signal  to  noise  ratio  spectrum  occur  for  lower  wavenvimbers . 
These  wavenumbers  are  the  spectral  components  that  contain 
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(arbitrary  unlta)  aa  a  function  of  wavenunber  (m“^)  for  cases 
of  a)hlgh  (2.90),  b)  Intermediate  (1.19)  and  c)  low  (.65) 
signal/noiaa  Tariance  ratio  (SNR).  Data  was  obtained  on 
January  19*  197b  at  White  Sands  Missile  Range  during  a  light 
snow  storm. 
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Figure  6:  Graphs  of  filtered  (lines)  and  unfiltered  [*) 
deviations  (Y'^'  (t^fRl)  from  equation  2.1.4)  in  the  mean 
aerosol  content  for  data  of  Fig.  3.  Note  the  scale  change 
for  large  SNR  values  (2.90). 
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velocity  inrorniilion.  Kip,.  i>  shows  a  decrease  in  the 
napnitude  and  bandwidth  of  the  filter  for  decrcuainp  SnR 
values.  Spectral  conponents  where  the  noise  contributions 
are  large  conpared  to  the  aerosol  inhosiogenei  ty  signal  are 
suppressed  an  order  of  magnitude  or  greater  by  this  process. 
Limitations  of  this  technique  are  evident  in  the  scatter  of 
the  estimated  signal  spectrum  to  noise  spectrum  ratios  for 
values  less  than  approximately  Ratios  less  than  this 

value  are  essentially  noise. 

Kig.  t>  shows  examples  of  filtered  and  unfiltered 
spatial  deviations  from  the  mean  aerosol  structure, 
calculated  from  Eq.  2.1.4.  Tne  filtered  aerosol  innomogen- 
eitiea  are  calculated  using  uue  uj^timal  transfer  functions 
shown  in  Fig.  5.  It  can  be  seen  from  this  figure  that  well 
defined  structures  corresponding  to  large  SNR  values  are  not 
substantially  changed,  but  random  fluctuations  are  consider¬ 
ably  reduced  by  this  filtering  procedure. 
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2.3  Velocity  Correction  for  Coherence  Decay 

A  nothod  is  developed  in  this  section  to  estimate  the 
mean  value  of  the  horizontal  wind  velocity  and  the  rms  wind 
speed  nodelinp  the  spatial  aerosol  inhonogencities  observed 
by  lidar  measurements  by  a  Gaussian  approximation.  A  deri¬ 
vation  of  the  functional  form  of  the  model  is  given  in 
Appendix  C.  The  final  form  of  this  model  includes  optimal 
filtering  characteristics  derived  in  the  previous  section  to 
maximize  the  information  in  the  profiles  of  aerosol  inhomo- 
geneities . 

The  spectral  velocity  component  perpendicular  to  the 
beam  pattern  u(k)  (see  Kig.  2)  is  calculated  from  the  coher¬ 
ence  of  lidar  measurements  at  various  times  and  spatial 
separations.  Random  changes  in  the  aerosol  inhomogeneity 
pattern  caused  by  turbulent  velocity  variations, decreases 
the  naxinun  value  of  the  coherence  with  increasing  time  lag 
and  lateral  separation.  Brigga  et  al.  (1950}*  Briggs  (196Ba, 
b)  and  Gosaard  (1969)  have  devised  first  order  approximations 
to  obtain  the  mean  ionospheric  drift  velocity  from  reflected 
radio  waves  in  the  presence  of  random  fluctuations.  They 
observed  that  the  correlations  and  cross  spectra  were  nearly 
Gaussian  in  time  and  spatial  separation. 

A  model  derived  in  similar  manner  to  that  of  Kunkel 
(1976)  is  used  to  estimate  the  mean  laterax  wind  opeea  of 
all  spectral  components  and  the  rms  wind  speed.  Appendix  C 
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preaents  a  derivation  of  the  fiinctional  form  of  the 
coherence  decay  with  increasing  temporal  and  spatial 
separations  for  an  ensemble  of  Gaussian  aerosol  inhomogen¬ 
eities  being  advected  by  a  Gaussian  velocity  distribution. 
The  functional  form  of  the  coherence  under  these  assumptions 
is  (Eq.  C.ll) 


2.3.1 


►/i 


where  k  represents  wavenumber  (m*^ ) 

u  is  the  mean  lateral  wind  velocity  component  (m/e) 
ax  is  the  average  lateral  spatial  separation  between 
Ildar  measurements  (m) 

is  the  ms  width  of  the  horizontal  extent  of  the 
spatial  aerosol  inhonogeneities  (m) 

C77is  the  standard  deviation  of  the  velocity 
distribution  (m/a) 

To  emphasize  the  spectral  components  where  the  signal 
spectrum  is  detectable  over  the  background  noise  and  to 
increase  the  statistical  significance,  a  weighted  average 
coherence  over  wavenumber  (Coh(AX,at))  is  computed.  The 
weighting  function  used  to  optimise  the  information  in  the 
signal  spectrum  over  the  noise  spectrum  is  the  square  of 
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the  optimal  tranal’cr  function  from  Kq.  since  factors 

of  the  optimal  transfer  function  occxir  in  spectral  estimates 
This  weighting  fur.ction  is  normalized  such  that  the  sum  of 
all  the  weights  is  unity. 

W(k)  =  K^(k)/  21  2.3.2 

♦•A. 


The  final  functional  form  of  the  model  is  the  weighted 
average  of  Eq.  2.3*1  over  all  wavenumbers 


where  Wik)  is  the  weighting  function  given  in  Eq.  2.3*2 
and  the  amplitude  factor  A  is  included  to  account  for 
finite  ratios  of  signal  to  noise. 

Experimentally  measured  weighted  coherence  calculations 
from  lidar  aerosol  inhomogeneity  profiles  are  fitted  to  the 
functional  form  of  Eq.  2.3.3  using  the  nonlinear  regression 
routines  provided  by  the  Madison  Academic  Computing  Center 
(Karquardt  1963).  The  four  parameters  regressed  against 


are 


A  the  amplitude  factor  to  account  for  finite  signal  to 
noise  ratios 

u  the  mean  lateral  wind  speed 

the  rns  horizontal  width  of  the  aerosol  inhomoge- 
nc i  t  ics 

C~j  the  standard  deviation  of  the  wind  speed  distribu¬ 
tion 

Experimental  data  points  for  nonlinear  regression  are 
chosen  about  the  maxiroun  weighted  coherence  value  for  each 
of  tne  three  lateral  separations  described  in  Section  2.1. 

An  example  of  these  weighted  coherence  measurements  for 
different  time  and  spatial  separation  along  with  a  least 
squares  fit  to  Eq.  2.3.3  is  shown  in  Pig.  7. 

An  initial  estimate  of  these  regression  parameters  ia 
needed  to  speed  convergence  of  the  nonlinear  regression  to 
the  optimal  values.  The  initial  estimate  of  the  amplitude 
factor  is  the  weighted  coherence  extrapolated  to  zero  time 
lag  for  Ax^O.  I 

A  -  CoF»l  0,  OJ  «  £  wu)  •  Aik)  2.3.U 

•  •4. 

The  weighted  coherence  are  assumed  to  have  a  maximum  value 
in  At  near  the  time  interval  needed  for  the  aerosol  inhomo- 
geneitiea  to  drift  the  lateral  separation  distanee  ^x. 
lEq.  2.1.5)*  This  assumption  is  valid  onlj  for  lateral 
velocities  much  greater  than  the  magnitude  of  the  wind 
speed  fluctuations.  Data  presented  in  this  thesis  sieet 
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Figur*  7t  Weighted  coherence  at  a  function  of  time  lag  for 
three  lateral  aeparationa  (^>0m(O  )•  Ax«36m(0  )) 

from  data  obtained  on  January  19*  197t)*  Linea  are  a  leaat 
aquarea  fit  to  Eq.  ?.3.3*  Relatively  poor  fit  forox-76B 
ia  a  reault  of  weighting  the  outaide  coherence  value  by  %. 
Theae  valuaa  are  weighted  beoauae  fewer  apeetra  are  uaed 
to  estimate  the  coherence. 
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this  requirement.  An  initial  estimate  for  u  under  these 

assumptions  is  the  apparent  lateral  velocity  u'  calculated 

from  the  time  laq  (^t  )  of  maximum  veiqhted  coherence  at 

m 

lateral  separation d x . 

u*  -  4kX/4t  «  2  R  sin(‘*/2)/^t  2.3.5 

m  in 

An  approximation  to  the  mean  square  width  of  the  Caussian 
aerosol  inhonogeneities  can  be  estimated  by  a  least  squares 
regression  of  weighted  coherence  values  for ^x»0  to  Eq.  2.3.3. 
The  rms  wind  speed  can  be  neglected  for  regression  of  data 
points  near4t*0.  Eq.  2.3.3  reduces  for'T^  approximately  zero 
to 

Coh(0,Jt)-  A  exp (-U * ^ A t^/2  C"  )  2.3.6 

where  is  the  only  remaining  unknown  initial  parameter  in 
this  equation.  The  initial  estimate  of  the  rms  wind  speed 
is  chosen  as  a  small  fraction  of  the  total  wind  speed  estimate. 

A  coarse  search  of  the  parameters  is  required  to  obtain 
a  close  estimate  for  a  true  minimum  and  not  a  local  minima. 

The  weighted  coherence  measurements  curves  are  regressed  to 
Eq,  2.3.3  for  each  parameter  while  holding  the  other  para¬ 
meters  constant.  This  procedure  is  iterated  until  the  change 
in  each  parameter  is  less  than  some  intermediate  error  toler¬ 
ance.  A  final  simultaneous  regression  on  all  parameters  is 
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perforned  to  establish  their  optir.al  values.  The  model 
parameters  to  be  compared  with  conventional  wind  velocity 
measurements  are  the  mean  lateral  velocity  component  u  and 
the  rms  wind  speed  cr^  . 

The  curve  fitting  procedure  described  above  is  used  to 
estimate  a  correction  from  the  apparent  velocity  caused  by 
turbulent  changes  in  the  aerosol  patterns.  -iTie  magnitude 
of  this  correction  depends  upon  the  magnitude  of  the 
maximum  time  lag  and  the  lateral  separation  distance.  This 
correction  can  be  as  much  as  10«20%  of  the  mean  lateral  wind 
speed  and  can  change  the  mean  wind  direction  estimaticn  3  to 
5  degrees. 

The  average  radial  velocity  (v^)  of  all  spectral 
components  is 


u 


1  v^(k) 

4.4  ^ 


2.3.7 


Radial  velocity  components  calculated  from  Eq.  2.1.9  are 
averaged  until  a  relative  phase  shift  in  the  cross  spectral 
estimate  greater  than  IbO  degrees  is  encountered.  This 
process  decreases  the  influence  on  the  mean  radial  wind 
speed  estimate  of  random  noise  phase  shift  in  higher  spectral 
components. 

Tne  horizontal  projection  of  the  mean  radial  velocity 


conponent  is 


V  =  Vp  cost  <9  )  2.3.0 

Data  presented  here  were  obtained  at  snail  elevation  angles 
(^10“)  and  therefore  cos(©>  "=  1.  The  nean  wind  speed 
(Vii^ar^  and  arinuth  direction  (  are  calculated  from 

the  following  formulas: 


V 

Uidar 


lu^  ♦ 


2.3.9 


^lidar  “  tan"\S/v)  2.3.10 

where  la  the  central  azimuth  angle  in  the  three  angle 
scan  of  Pig.  2. 
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2.4  Lidar  Uncertainty  Estimates  in  Wind  Measurements 

An  estimate  in  the  experimental  uncertainties  is  needed 
to  Judge  the  reliability  of  any  measurement .  /knowledge  of 
the  factors  contributing  to  measurement  uncertainties  can  be 
used  to  improve  experimental  procedures.  Lidar  wind  measure¬ 
ment  uncertainties  are  based  upon  the  discreteness  of  lidar 
data.  Lidar  measurements  of  aerosol  profiles  are  digitized 
for  ease  in  data  manipulation  and  represent  measurements 
averaged  over  a  range  interval  of  l^n  (see  Section  3>1}* 
i*hese  measurements  are  acquired  at  discrete  time  intervals 
limited  by  the  laser  pulse  repetition  rate. 

The  wind  speed  and  azimuth  angle  ( ^ )  between  the  wind 
azimuth  )  and  the  average  lidar  azimuth  (  (i.e., 

,c*n.  ^  can  be  estimated  from  the  approximate 
formulas  for  the  velocity  components  u  s  ^x/et^  and  v  s 
Ay/^it^  by 


^  !  Wu*«v  ^ 


2.4.1 

2.4.2 


where  ax  is  the  average  lateral  separation 

Ay  Is  the  rsdial  spatial  change  (units  of  15m)  of 


maximum  correlation 
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At  is  the  time  lag  for  maximum  correlation  which  is 
n 

the  approximate  time  interval  for  the  aerosol 
inhomogeneities  to  drift  the  lateral  distance  a*. 
i;^uanti zation  uncertainties  in  lidar  wind  velocity 
measurements  for  fixed  value  of  a x  can  be  estimatea  f  ru«*« 


following  equations.  The  maximum  azimuth  quantization 
uncertainty  estimated  from  Eq.  2.4.2  for  changes  in 


2.4.3 


The  maximum  quantization  uncertainty  in  the  wind  speed 

is  estimated  from  Eq.  2.4.1  for  discrete  changes  in  the  reaxi> 

mum  time  lag  ^t.. 

n 

S  V  .■  .l.)j  J.4.1, 

A  schematic  diagram  of  quantization  uncertainties  in 
speed  and  azimuth  is  shown  in  Pig.  b.  The  uncertainties  are 
plotted  as  functions  of  wind  speed  and  azimuth  angle  (0^) 
with  respect  to  the  average  lidar  azimuth  direction.  The 
graph  represents  lidar  quantization  uncertainties  fori9x^7.9B 
and  a  1  hs  laser  pulse  repetition  rate. 

Uncertainties  in  the  lidar  rms  wind  speed  measurements 
are  estimated  by  the  nonlinear  regression  (see  section  2.3). 
The  nonlinear  regression  routines  estimates  the  9S%  confidence 


limits  of  the  regrcsaion  parameter  Cj  .  The  magnitude  of 
these  values  are  used  for  lidar  nas  wind  speed  uncertainties 


LV. 

I  ?  i  S  I  .  A  b  to  ’ 


SPFEO  (M/S) 


Figure  6i  Maximun  lidar  quantization  uncertainties  (Eqs.  2.U.3 
and  2.U.U)  for  ^x-U7.9ni  and  a  1.0  Hz  laser  pulse  repetition 
rate.  Solid  lines  are  Piaxinure  speed  quantization  uncertainties 
(m/ase)  and  dashed  lines  are  maximum  azimuth  direction 
uncertainties  (degrees!. 
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3.1  The  Wisconsin  Lidar  System 

The  Wisconsin  Lidar  system  is  operated  in  a  monostatic 
configuration.  A  ruby  laser  pulse  is  transmitted  into  the 
atmosphere  and  the  energy  output  is  measured.  The  laser 
radiation  scattered  through  IbO  degrees  by  aerosols  and 
molecules  is  detected  by  a  receiver  telescope  aligned  parallel 
to  the  laser  axis.  A  IP  nm  bandpass  interference  filter  is 
used  in  the  receiver  optics  to  suppress  the  background 
radiation. 

The  field  of  view  of  the  telescope  is  variable  and  full 
overlap  of  the  receiver  field  of  view  and  the  laser  beam 
divergence  occurs  at  approximately  1  km.  A  red  sensitive 
photomultiplier  converts  the  backscattered  laser  photons  into 
an  electrical  signal.  The  photomultiplier  signal  is  loga¬ 
rithmically  amplified  to  compress  the  dynamic  range  and 
digitized  at  a  10^5i^  rate  by  a  Eiomation  model  1010  10-bit 
analog  to  digital  converter.  The  lOKhz  rate  yields  a  range 
resolution  of  15  m.  The  Wisconsin  lidar  system  parameters 
are  listed  in  Table  1.  A  block  diagram  of  the  system  (Kig.  9) 
is  also  included. 

A  PDF  11/40  minicomputer  is  used  for  data  preprocessing 
and  displaying.  Individual  lidar  returns  are  corrected  for 
the  range  square  attenuation  by  adding  twice  the  natural 
logarithm  of  the  range  to  each  data  point.  Returns  are 
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Figure  9t  Lidar  ayataa  block  diagraa. 
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norr.alired  by  aubtructinp  the  natural  lofarit^js  of  the 
output  energy  for  each  laser  pulse.  Real  tir.e  corrected 
lidar  returns  ‘•'I*  «*re  then  displayed 

on  a  CRT  monitor. 

The  laser  telescope  system  is  mounted  on  a  motorized 
base  which  allows  the  positioning  of  the  lidar  in  azimuth 
and  elevation  by  the  PDF  11/UO  computer.  Angle  readout 
encoders  positioned  on  the  base  detemine  the  angular 
coordinates.  These  angles  are  recorded  on  magnetic  tape 
along  with  the  corrected  lidar  returns,  time,  date  and 
laser  pulse  energy. 

The  Wisconsin  lidar  system  can  be  operated  from  an 
electronics  trailer  for  field  experimentation.  The  trailer  is 
equipped  with  power  hookup,  air  conditioning  and  a  cooling 
system  for  the  laser.  The  computer,  laser  power  supply  and 
lidar  base  are  equipped  with  shock  mounts  to  minimize  vibra* 
ti  CDS  in  transit.  The  lidar  base  is  extended  onto  a  lowered 
tailgate  during  measurements. 


3.2.  Tower  Wind  ^feasurcncnt  Inetra":entflti  on 


Independent  wind  measurcrienta  in  exporir.onta  conducted 
at  !’.adlson,  Wisconsin  were  obtained  for  co/nparlson  with 
Ildar  derived  quantities,  k  Gill  four  blade  anemometer 
IR.rt.  Young  model  P1002)  was  mounted  at  a  height  of  77n 
from  the  base  of  a  radio  tower  located  2,7  lui  west  of  the 
Keteorologj  and  Space  Science  Building  at  the  University 
of  Wisconsin,  Madison  (see  Fig.  10).  A  12«bit  A  to  D 
converter  is  used  to  dlg*ti7.e  the  wind  speed,  azimuth  and 
elevation  angles  of  the  bivane  sequentially,  each  at  a 
1.0  Hz  rate.  Data  transfer  to  the  PDF  11/UO  based  data 
logging  system  (Section  3.1)  is  done  via  a  telephone  modem 
link.  Tower  wind  data  are  stored  on  digital  magnetic  tape 
along  with  concurrent  Ildar  measurements. 

A  number  of  factors  contribute  to  the  uncertainty  in 
comparison  of  tower  anenoneter  wind  measurements  wtti  lidar 
wind  measurements.  Sethuhanan  and  Brown  (1976)  estimate  that 
prop  response  errors  in  wind  speed  measurements  are  less  than 
approximately  2%.  Aerodynamic  lifting  of  the  bivane  tail 
produces  errors  in  the  angular  positioning  smaller  than  two 
degrees,  cf.  Pendergsst  (1975).  The  supportive  tower 
structure  can  also  influence  wind  measurements.  Wind  speed 
measurements  can  deviate  up  to  from  the  mean  speed  for 
sensors  downwind  from  the  tower  sccording  to  Ituni  and  Barad 


(1S70).  Upwind  sensors  can  te  alTected  up  to  b%  in  a 
similar  r.onner  according  to  Angell  and  Bernstein  (1976), 
Tower  structure  influences  upon  unenomter  wind  r.easurer.er.ts 
are  estimated  to  be  much  smaller  than  these  extreme  values 
in  this  study.  Lidar  and  tower  anemoroter  wind  measurements 
represent  averages  over  different  sample  volumes  and  there¬ 
fore  different  statistical  sampling.  Tower  wind  measurement 
uncertainties  are  based  upon  the  statistical  fluctuations 
in  these  measurements. 

Estimation  of  tower  anemometer  wind  velocity  uncertain¬ 
ties  depend  upon  the  correlation  of  these  measurements  in 
time.  Error  analysis  based  upon  a  statistical  description 
is  detailed  in  Appendix  D.2.  An  independent  measurement  of 
the  average  wind  field  is  assumed  to  be  made  during  a  time 
period  equal  to  the  Eulerian  time  scale  (  Z*  )  for  stationary 
isotropic  turbulence.  The  uncertainty  in  the  measure  of  the 
mean  value  of  a  wind  parameter  is  estimated  by  Eq.  D.2  .  3 
to  be 

C-/  :  O; /(«,/'. 

where  is  the  standard  deviation  of  the  parent  population 
for  parameter  p 

ia  the  uncertainty  in  the  mean  value  of  the 
parameter  p 


45 


3.3  V.’hite  Snnds  Kiasile  hanr.c  inatruncntstion 

Intlcpetident  wind  irenaurcncnla  for  coi'.purison  to  Ildar 
wind  ncaaurcnents  were  obtained  for  cjcporir.cnta  conducted  at 
the  White  Sands  Missile  Range,  New  Mexico.  A  map  of  the 
experimental  site  at  White  Sands  Launch  Complex  |r3b(LC3(>)  is 
shown  in  Kig.  li.  Wind  velocity  measurements  from  an  instru¬ 
mented  meteorological  tower  and  radar  tracked  pilot  balloons 
are  used  in  this  thesis  for  comparison.  Experiments  were 
conducted  from  December  7-17,  1*^77  and  January  7-27,  1V70. 

The  meteorological  tower  at  the  LC36  sice  was  instru¬ 
mented  for  wind  measurements  witn  the  following  equipment! 
cup  anemomters  were  at  eight  height  levels  of  the  tower 
i7.b,  23,  30,  b3»  b9,  91,  122,  152  m).  Measurements  from 
the  cup  anemometers  were  averaged  every  15  seconds  and  re¬ 
corded  on  paper  tape.  Two  U-V-W  anemometers  (R.H.  Young 
model  27U03)  were  located  at  heights  of  lU  and  137  n.  Wind 
measurements  from  these  two  anemometers  were  acquired  by  the 
same  data  transfer  system  described  in  Section  3.2  and  stored 
along  with  lidar  measurements. 

Wind  measurements  from  automatically  tracked  pilot  balloons 
were  used  for  comparisons  with  Ildar  wind  profile  measurements 
above  the  tower.  Pilot  balloon  measurements  are  vertically 
averaged  over  heights  of  100  m.  Averaging  pilot  balloon 
measiirements  in  this  manner  yields  the  same  vertical  averaging 


Figure  lit  Map  of  Launch  Complex  #36(LC36),  axparimental 
alia  at  White  Sands  Hiasile  Range,  New  Mexico.  Lidar  and 
toner  positions  are  marked  L  and  T. 
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as  lidar  wind  measurements  and  smooths  over  variations  in 
wind  velocities  due  to  nutations  in  tne  radar  dish.  Pilot 
balloon  rr.easurcnonts  were  obtained  every  10  minutes  ior  the 
data  presented  here. 

Uncertainties  in  average  pilot  balloons  wind  measurements 
can  be  estimated  from  inaccuracies  in  the  average  spatial 
displacement  of  the  balloons  in  time.  Displacement  errors 
cati  uccur  beca..se  wx  erratic  motion  of  these  balloons  due  to 


short  time  scale  turbulent  gusts  and  aerodynamic  forces  upon 
the  balloon  during  ascent,  cf.  Rider  and  Amendariz  (196ti), 
Rogers  and  Camnitz  (19&6).  Displacement  errors  also  occur 
because  of  positioning  inaccuracies  due  to  equipment  resolu* 
tion,  see  Schaefer  and  Doswell  (197h).  Pilot  balloon  wind 
velocity  uncertainties  are  estimated  in  this  thesis  from 
positioning  errors  in  space. 

Pilot  balloon  wind  velocity  measurement  uncertainties 
are  derived  in  Appendix  D.l  for  spatial  positioning  errors. 
The  maximum  uncertainty  in  each  velocity  exponent  ^  , 


^  if  ^ 

U  —TZ - 


3.3.1 


where  R  is  the  average  slant  range  of  the  pilot  balloon  during 
the  measvirement 

(  R  is  the  uncertainty  in  the  slant  range  measurement 
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U.O  Experimental  Kesults 

Experiments  were  conducted  in  Kadison,  Wisconsin  and 
at  the  ’..“hite  Sands  Missile  Ranpe,  New  Mexico  to  test  the 
theory  presented  in  Chapter  2.  Aerosol  density  inhomoge- 

neit>  scale  si:es  from  60  m  to  1  km  were  found  to  be  efficient 
wind  tracers.  These  experimental  results  represent  the  mean 
motion  of  the  aerosol  inhomogeneities  over  a  five  minute 
period.  Uncertainty  estimates  of  the  wind  velocity  are  listed 
for  both  lidar  and  comparisons  measurements.  Results  show 
that  lidar  can  measure  the  spatial  and  temporal  variations 
in  the  wind  velocity  field. 


U.l  liadison  Experiments 

Conparisons  of  lidar  and  tower  anemometer  wind  velocity 
measurements  were  conducted  in  Kadison,  Wisconsin  from  April 
through  August  1977*  Remote  lidar  measurements  were  obtain¬ 
ed  over  the  western  section  of  Kadison  (see  Fig.  10).  The 
Wisconsin  lidar  system  was  located  on  the  ninth  floor  of  the 
University  of  Wisconsin  Meteorology  and  Space  Science  Building. 
Results  of  these  experiments  are  for  convective  days  with 
well  defined  inhomogenoities  in  the  mean  aerosol  structure. 
Aerosol  inhomogeneities  are  of  sufficient  magnitude  under 
these  conditions  to  be  detectable  above  the  background  noise 
of  the  system.  A  signal  to  noise  ratio  (SM))  defined  from 
Eq.  2.2.5  of  greater  them  0.5  occurred  for  the  data  present¬ 
ed  in  this  section. 

Lidar  profiles  of  aerosol  inhomogeneities  were  obtained 
in  the  three  angle  azimuth  scan  described  in  Section  2.1. 

Lidar  measurements  were  acquired  with  a  1.0  degree  azimuth 
beam  separation  at  an  elevation  angle  of  2.3  degrees  for 
each  day.  The  tower  anemometer  location  (shown  in  Fig.  10) 
was  located  near  the  center  of  the  azimuth  scan  pattern.  A 
data  segment  length  of  960  m  (64  data  points)  centered  about 
the  tower  distance  (2.7  km)  was  used  in  analysis  of  wind 
information.  The  average  lateral  beam  separation  (Eq.  2.1.5) 
was  47>9  m.  The  period  between  successive  lidar  aerosol 
inhonogeneity  profiles  was  limited  by  laser  cooling  requirements 
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to  1.1  seconds.  Lidar  wind  measurements  represent  five 
minute  time  averaces  over  a  sample  volume  960  m  long  by 
47.9  m  wide  by  17  m  high. 

Lidar  aerosol  measurements  were  analyzea  xor  wind 
velocity  information  by  the  method  outlined  in  Chapter  2 
and  compared  with  concurrent  tower  anemometer  wind  measure¬ 
ments.  The  experimental  time  duration  of  each  day  varied 
between  20  to  45  minutes  so  that  a  number  of  five  minute 
wind  velocity  averages  could  be  computed.  Uncertainties  in 
lidar  wind  raeasuroraents  are  estimated  according  to  Section 
2.4. 

Uncertainties  in  tower  anemometer  measurements  (see 
Section  3.2)  are  estimated  as  follows.  The  Eulerian  speed 
autocorrelation  is  calculated  from  the  tower  measurements 
over  the  entire  time  duration  of  the  experiment.  The 
Eulerian  tine  scale  is  estimated  as  the  tine  integration  of 
the  speed  autocorrelation  function  to  a  tine  where  this 
value  is  approximately  zero.  Uncertainties  in  the  measure¬ 
ment  of  the  average  tower  speed  and  azimuth  direction  f  (1. 
are  estimated  by,  cf.  Eq.  3.2.1. 

5  >  CTf  '  ) 

.  »  -  i'L 

‘  ^•/“T  ) 


4.1.2 
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where  2^  is  the  Kulcrian  tine  scale 

T  is  the  tine  duration  of  the  entire  experiment 

s  the  tower  nns  wind  speed  calculated  over  the 
entire  experimental  tine  duration 
Oj  ia  the  uncertainty  in  the  average  tower 

wind  speed  measurement  V 
is  the  rms  tower  wind  azimuth  direction 
calculated  over  the  entire  experimental  time 
duration 

.  C^-  is  the  uncertainty  in  the  average  tower 

wind  azimuth  direction 

The  ms  wind  parameters  determined  for  the  experimental 
duration  (T)  are  assumed  to  approximate  the  standard 
deviation  of  the  parent  populations.  Uncertainties  in  the 
measurement  of  the  five  minute  averaged  tower  rms  wind  speed 
estimates  are  calculated  as  the  standard  deviation  of  these 
measurements  from  the  ms  wind  speed  calculated  for  the  entire 
time  period  ( 


U.1.3 


where  is  the  number  of  five  minute  average  wind 
measurements  for  each  experiment 
is  the  tower  rms  wind  speed  measured  over  a  five 
minute  average 
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is  uncertainty  in  the  tower  eatination  of 
for  each  experinent 

These  uncertainties  in  tower  ancr.or.eter  ineasurer.cntc  are 
listed  in  Table  2. 

The  results  of  these  experinenta  are  stunnarised  as 
time  histories  of  wind  speed,  azimuth  direction  and  ma 
wind  speed  of  lidor  and  tcver  anemometer  measurements  in 
Figs.  12-17.  Wind  speeds  and  directions  are  calculated  by 
Eqs.  2.3.9  and  2.3»30.  The  average  difference  between  lidLar 
and  Luwer  mc a sui  w*  whe  .MCttn  wind  speed  and  azimuth 

direction  as  well  as  the  3ta:»dnrd  deviations  of  these  errors 
are  listed  in  Table  2  for  each  experimental  day,  iTiese  are 
calculated  from  the  following  formulas: 
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Pigur*  12t  Conpariaons  of  lidar  (♦«  aolid)  and  towar  (0, 
daahad)  a)  apaad,  b)axliButh  (froia  north)  and  e)  ^ 

▼eraua  ti«a  for  22  April  1977.  Data  obtainad  In  Kaoiaon. 
All  maaaurananta  made  at  the  tovar  anemomatar  haight  (77m). 
RMS  wind  apaad  for  the  second  data  point  is  not  plotted 
because  the  coherence  maxima  wars  nearly  as  low  as  those 
correaponding  to  a  5S  random  noise  coherence  (Eq.  A2.2). 
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Tine  I  COT  I 


Figure  17 1  Similar  to  Figure  12  for  26  July  1977.  Non-linear 
regreaelon  eonTerged  to  lover  bound  In  eigma  due  to  lov  SNR 
for  aeeond  to  laat  point.  Thla  data  point  la  not  plotted. 

(aee  page  vi  ). 
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where  those  parainotcrs  are  described  in  Table  2.  ‘i'he 
ovoracc  lidar  quantization  Ui<certaintios  ore  also  listed  in 
Table  2. 

Klgs.  lb,  19,  and  20  are  conparisons  of  lidar  and  tower 
anenometer  measurenonts  of  the  wind  speed,  azimuth  direction 
and  the  ms  wind  speed,  respectively.  Average  uncertainties 
in  both  lidar  and  tower  measurements  for  each  experiment  are 
also  shown.  These  graphs  show  that  lidar  wind  measurements 
are  consistent  with  concurrent  tower  wind  measurements. 

Some  limitations  of  this  procedure  for  remote  lidar  wind 
measurements  were  discovered  while  investigating  measurements 
with  large  discrepancies.  The  process  of  removing  the  mean 
aerosol  structure  (Eqs.  2.1. 3*2.1. 4)  bias  against  detecting 
small  radial  displacements.  Aerosol  inhomogeneities  do  not 
have  time  to  drift  the  radial  distance  between  data  points 
il5  m)  in  the  time  for  them  to  drift  laterally  between  beams. 
This  occurs  for  radial  wind  speeds  typically  less  than  a 
meter  per  second.  The  large  discrepancies  in  azimuth  on 
Hay  23  and  July  2b  can  be  explained  by  the  overestimation  of 
lidar  radial  wind  speeds.  Lidar  azimuth  quantization 
uncertainties  are  large  in  these  instances. 

Extremely  low  values  of  lidar  measured  rma  wind  speeds 
on  April  22,  June  23  and  July  27  were  found  to  have  low  SNR 
values  I'he  coherence  maxlAS  in  these  cases  were  the 

sane  order  of  magnitude  as  the  3b  probability  coherence  for 


LIDflR  SPEEDS  (H/SEC) 
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TOWER  SPEEDS  (M/SEC ) 


figur*  Ibt  Comparlaoni  of  Ildar  and  averaga  tower  anemoneter 
wind  apeeds  for  April  22  (□  ),  May  23(0)>  June  20  (^), 

Juno  23  (♦),  June  21  (0)»  And  July  27  (♦),  1977.  Data 
aognanta  warn  obtained  on  elaarf  convaetiTa  days.  Lidar  and 
tower  uneartaintias  are  ahown  for  each  axpariaiental  day. 


UJ  40 
Ql 


lui 


10 


v 


13*  I  *0 


ai 

a: 

Q  340 


TOWER  nZIMUTH  (DEGREES) 


Figure  19 1  Slnilar  plot  to  Figure  18  for  the  etimuth  angle 
(degrees  from  North).  Lidar  and  tower  azimuth  uneertaintiea 
are  alao  ahown.  Kotei  Data  on  May  23>  June  23*  ^nd  June  27 
are  plotted  for  the  oppoaite  direction,  i.e.,  160  degree 
difference . 
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Pigur*  20 t  Slnllar  plot  of  rma  wind  speed  for  the  date  of 
Figure  l8«  Lidar  and  tower  uncertainties  for  each  ezperinental 
*^•7  slso  shown.  Valuea  where  the  non-linear  regression 
eonvergeu  i,o  vary  lo*  Tslues  because  of  low  SKP  values  (see 
pageOi)  are  not  plotted. 
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random  fluctuations  predicted  by  Eq.  i'nese  coherence 

maxima  -..ere  found  not  to  decrease  with  increasing  lateral 
separation  as  in  Eig.  7.  V.'ind  speeds  and  directions  calcu¬ 
lated  from  this  data  were  found  to  be  consistent  with 
previous  measurements,  however  the  nonlinear  regression 
converged  to  very  low  values  to  minimize  the  residue  in  the 
regression.  These  measurements  of  the  ms  wind  speed  are  not 
plotted  because  of  the  large  uncertainty  in  the  regression  of 
this  parameter. 


U.2.  White  Snnds  Kxi’criments 

A  field  experir.ent  wbs  conducted  in  conjunction  with 
the  Atmospheric  Science  Laboratory  at  the  V'nlte  Sands 
Kissile  Range,  New  Mexico,  henote  lidar  wind  measurement 
experiments  were  conducted  on  the  dates  December  7-17,  1977 
and  January  7-22,  197tJ.  Fig.  11  is  a  map  of  the  experi¬ 
mental  area.  The  Wisconsin  lidar  system  was  located  2.7  km 
south  of  the  Launch  Complex  #3^  (LC36)  meteorological  tower. 
Instra-itntation  for  tower  wind  measurements  are  described  in 
Section  3.3* 

The  experimental  procedure  for  obtaining  remote  lidar 
wind  measurements  was  similar  to  that  used  in  the  Madison 
experiments  (Section  U.l).  The  three  angle  aiimuth  scan 
described  in  Section  2.1  is  used  with  an  azimuth  separation 
between  lidar  return  aerosol  profiles  of  one  degree.  The 
lidar  elevation  angle  was  positioned  at  6  degrees  for  height 
profiles  of  wind  velocity.  The  length  of  the  data  segments 
(OU  points,  960  m),  laser  firing  period  (1.1  seconds)  and 
time  duration  for  each  measurement  (5  minutes)  are  the  same 
as  those  used  in  Section  U.l. 

Height  profiles  of  wind  velocity  correspond  to  lidar 
measurements  at  different  ranges  along  the  propagation  path. 
Average  height  levels  of  lidar  measurements  are  calculated 
froei  Eq.  2.1.6.  The  volume  over  which  lidar  wind  measurements 
ere  averaged  varies  because  the  lateral  separation  increases 
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with  Increasing  ranee.  The  average  lateral  separation  varies 
from  21  n  at  a  slant  range  of  1.2  km  to  99  m  at  5*7  km.  The 
vertical  and  lonciiudinal  lengths  of  the  sampling  volume  are 
100  m  and  9b0  m,  respectively,  for  all  ranges. 

Wind  velocity  measurements  from  lidar  data  are  calculated 
using  the  theory  in  Chapter  2.  Wind  velocity  height  profiles 
are  measured  using  nonoverlapping  data  segments.  The  lowest 
height  level  is  limited  by  the  range  for  complete  overlap 
between  the  laser  beam  divergence  and  the  field  of  view  of  the 
receiver.  Uncertainties  in  lidar  wind  measurements  are  esti* 
mated  according  to  Section  2.U. 

Measurements  presented  in  this  section  were  obtained  from 
1620-1720  MST  on  January  19,  1970.  Very  light  snow  with  large 
scattering  cross  section  in  the  visible  was  observed  to  fall 
throughout  the  period  (visibility  0  kn).  Large  variations  in 
the  return  lidar  signal  were  observed  to  a  range  of  6  km. 

Lidar  observations  at  higher  elevation  angles  showed  multiple 
cloud  levels  with  the  lowest  at  approximately  bOO  m.  This  data 
was  chosen  to  illustrate  the  capabilities  of  lidar  in  cases  of 
high  signal  to  noise  ratios  (SKR). 

Results  from  this  data  are  shown  in  Figa.  21-26.  Each 
graph  consists  of  height  profiles  of  wind  speed,  azimuth 
direction  and  rms  wind  speed  for  lidar  and  pilot  balloon 
measurements.  Pilot  balloon  measurements  are  plotted  for 
balloons  before,  during,  and  after  lidar  measurements,  when 


16U5-1<>50HST.  Pilot  balloons  are  at  1633-1637  (□  )»  16U3-16U7 
and  1655-1637  {A  )  for  January  19*  1976  data  at  White  Sande. 
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applicable.  These  balloon  measurements  are  averaged  for 
comparisons  with  lHar  wind  profiles. 

Keasurementa  of  the  mean  wind  speed,  direction,  and 
rns  speed  from  the  U-V«W  anemometer  at  the  137  m  tower 
level  at  LC3t>  are  also  plotted  in  Figs.  21«26.  The  north 
south  component  was  found  to  be  consistently  smaller  than 
lidar  and  pilot  balloon  measurements.  Comparison  of  this 
component  with  the  cup  anemometer  at  the  1^2  m  level  indi¬ 
cates  that  the  reliability  of  the  wind  sensor  in  this 
direction  is  questionable. 

Differences  between  lidar  and  average  pilot  balloon 
wind  measurements  are  listed  in  Table  3  along  with  the 
mean  and  rms  speed  and  azimuth  errora  for  each  height  level. 
Average  lidar  quantization  and  pilot  balloon  uncertainties 
for  each  level  are  also  known. 

Pigs.  27*32  ere  time  histories  of  lidar  and  pilot 
balloon  mean  wind  speed  and  azimuth  measurements  for  each 
height  level.  Pig.  27  also  shows  the  wind  speed  and  direction 
measured  by  the  U-V-W  tower  anemometer  at  the  137  m  level. 

The  large  variations  in  the  wind  direction  are  explained  by 
the  instrumental  difficulties  in  the  north  south  wind  speed 
sensor.  Broad  coherence  peaks  with  multiple  maxima  were 
observed  for  the  59U  ■  lidar  wind  measurement  at  lb3b  MST. 

The  maximum  value  occurred  for  a  time  lag  smaller  than 
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Figur*  21 t  Tims  history  of  lidsr  {*),  snd  pilot  bslloon  (0  ) 
nsasursmsnts  st  ths  127b  Isssl  for  data  obtainsd  on  January  19, 
1976  at  Whits  Sands  Mi sails  Rangs.  Tovar  ansmomstsr  nsaaursmsnts 
(*)  St  ths  137b  IsTsl  srs  also  shown.  Ths  Isrgs  Taristion  in  ths 
towsr  asiButh  is  eaussd  by  diserspaneisa  in  ths  north  south 
eomponsnt,  sss  psgs  7</  . 
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TIME  (MST) 


Figiire  29 1  Sinilar  graph  for  the  320a  level 
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Figure  32 1  Siailar  graph  of  Figure  27  for  the  591ui  lerel. 
The  firat  data  point  exhibited  broad  eoherenee  peaks  with 
Multiple  Maxima.  High  wind  speed  corresponded  to  a  shorter 
time  lag  than  expected,  aee  page  • 
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expected  fran  pilot  balloon  measurements,  producing  higher 
lidar  wind  speed.  Vertical  displacement  of  snowfall  varia* 
tions  could  be  significant  for  the  large  lateral  separation 
(99  m)  at  this  height. 

Figs.  33  a^d  3U  ^re  comparison  of  mean  wind  speed  and 
atinuth  direction  between  lidar  and  pilot  balloon  measurements. 
The  minimum  and  maximiua  lidar  quantisation  and  pilot  balloon 
uncertainties  are  also  shown.  These  limiting  values  are 
graphed  since  data  points  at  different  height  levels  are 
within  close  proximitj  of  each  other. 
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PIBflL  SPEED  (M/SEC) 


Figur*  33 t  Coaparlsona  of  Ildar  and  averaga  pilot  balloona 
apoad  (m/a)  for  aaeh  haight  laral  oomputad  (□  -127m,  O  •227m, 
d -320b,  ♦  -410b,  ^  501b,  If  594*)  •  Tbaaa  rasulta  ara  froB  data 
obtainad  on  January  19*  1978  at  Uhlta  Sanda  Kiaalla  Ranga. 
KlniauB  and  naxlBUB  Ildar  quantisation  and  pilot  balloon 
uneartalntlas  froB  Tabla  3  nra  plottad  alnea  Baasuraaanta  ara 
In  elosa  proslBlty  with  aaeh  othar. 


120  ISO 

PIBflL  AZIMUTH  (DEGREES) 


Figure  34t  Sinllar  plot  to  Figure  32  for  the  azimuth  direetiona 
of  ell  levela  (sane  symbols)  for  January  19  data.  Minimum  and 
maximum  uneertainties  fr<Mi  Table  3  are  also  shown. 


5.0  Conclusions 

A  procedure  has  been  developed  for  remote  measurements 
of  the  wind  velocity.  The  motion  of  naturally  occurring 
aerosol  density  inhoraogeneitles  are  measured  by  Ildar  to  obtain 
average  wind  speed,  direction  and  rma  speed  estimates.  The 
Past  Fourier  Transform  (PPT)  has  been  Implemented  In  this 
procedure  to  Increase  computational  efficiency. 

An  Improvement  In  spatial  filtering  of  Ildar  data  has 
also  been  developed.  The  ratio  of  the  signal  to  noise 
spectral  characteristics  are  used  to  design  an  optimal  linear 
filter  based  upon  the  Ildar  measurements  themselves.  Spectral 
components  irttere  the  background  noise  is  large  compared  to  the 
aerosol  Inhomogeneity  signal  are  suppressed  one  to  two  orders 
of  magnitude  with  respect  to  those  with  large  signal  eomponents 
(see  Fig.  5).  Aerosol  Inhomogeneity  scale  sixes  seen  from  this 
flgxire  range  from  bO  m  to  1  km.  These  scale  sixes  carry  wind 
Information  detectable  with  the  present  system.  The  ratio  of 
the  signal  spectrum  to  the  noise  spectrum  (Pig.  4)  show  high 
variability  between  successive  spectrsl  eomponents  for  values 
ssMller  than  approximately  0.5«  Lidar  aerosol  inhomogeneity 
profiles  with  SUB  less  than  0.5  are  therefore  dominated  by  the 
background  noise  of  the  present  system.  Wind  information  in 
these  instances  are  not  recoverable. 

Estimation  of  the  mean  wind  velocity  in  the  turbulent 
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planetary  boundary  layer  has  been  improved  by  including  a 
correction  for  the  spatial  and  temporal  coherence  decay, 
hesults  from  lidar  wind  measurenents  follow  the  temporal  and 
spatial  characteristics  of  the  boundary  layer  flow.  Comparisons 
of  lidar  wind  measurements  with  independent  wind  measurements 
shown  in  Kigs.  Id,  19*  33*  and  34  are  generally  within  the 
estimated  experimental  uncertainties. 

A  measiu'e  of  the  turbulent  dispersal  of  the  aerosol 
inhomogeneities  can  be  estimated  from  the  coherence  decay  with 
increasing  spatial  and  temporal  separation.  The  model  described 
in  Section  2.3  ia  used  to  estimate  the  turbulent  rms  wind 
speed  which  cause  the  aerosol  inhomogeneities  to  change  in 
time.  Kunkel  (197d)  has  used  similar  measxirements  to  estimate 
the  turbulent  energy  dissipation  rate  for  all  scales. 

Lidar  rms  wind  speed  measurements  are  shown  in  Kig.  20 
to  be  consistent  with  concurrent  tower  anemometer  estimates. 

A  measure  of  the  reliability  of  lidar  wind  measurements 
can  be  obtained  from  the  ratio  of  the  measured  error  between 
lidar  and  independent  measturements  to  the  uncertainty  in  these 
measurements.  This  is  expressed  by  the  following  formulas. 

5.1 

5.2 
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wnere  41/  is  the  wind  speed  difference  between  the  lidsr  and 
independent  wind  measurements  defined  in  Eq. 
d/  is  the  azimuth  direction  difference  between  iiuar 
and  Independent  wind  measurements  defined  in  £q. 

is  the  rms  value  of  the  lidar 
and  independent  wind  direction  uncertainties. 

.  *  ia  the  rms  value  of  the  lidar 

0  * 

and  independent  wind  speed  uncertainties. 

is  the  ratio  of  the  magnitude  of  the  speed  difference 

to  the  rms  speed  uncertainties. 

6^  is  the  ratio  of  the  magnitude  of  the  azimuth 
difference  to  the  rms  azimuth  uncertainties. 

Graphs  of  these  ratios  (  )  are  shown  in  Fig.  33 

as  a  function  of  the  SNR  estimate  (Eq.  2, 2, 5),  The  average 
values  of  and  along  with  the  standard  deviations 
are  shown  in  the  intervals  SNR  £  1.5*  1.3  ^  SNR  £  2.3,  SNR ^ 


2.3. 


Values  larger  than  one  standard  deviation  from  the  mean 


represent  a  statistical  probability  (33%)  that  the  error  is 


larger  than  the  uncertainty.  Thit  fig\u*e  shows  the  average  value 
of  and  aa  well  as  the  probability  for  larger  errors  than 
the  uneertaintiea  decreases  with  increasing  SNR  values. 


Large  discrepancies  of  some  data  points  can  be  explained 
by  the  limitations  of  the  method  and  equipment.  The  circled 


ee 


t 


Figur*  35t  Orapha  ot  ,&j  (•••  Eqa.  5*1  and  5.2)  vtaieh  ara 
funetioDi  of  tha  SXR  (Eq.  2.2.5)  ealoulatad  froai  Kadiaon  and 
Uhita  Sanda  data.  Sxparlnantal  daya  ara  April  22  (O  ),  May  23  (0), 
Juna  20  (4  )»  Juna  23  (♦)•  Juna  27  ( X )  and  July  26  (Q),  1977. 

.  Whita  Sanda  data  pointa  ara  plottad  for  tha  halght  laaala 

(  -  127*  (t  )^  227*  (X)»  320n  (X)»  410*  (-di),  501*  (*  )  and 

r  :  59iMB  (9).  Cirelad  and  boxad  data  pointa  ara  daaeribad  in  tha 

>  j  taxt  on  paga  97  . 


data  points  in  Fig.  represent  r.easurenents  where  the 
radial  wind  speeds  are  snail  (^In/s).  Large  azimuth 
differences  between  lidar  and  tower  neasurenents  are  due  to 
larger  lidar  radial  velocity  estimates  (see  discussion  on 
page^f).  The  boxed  data  points  in  Pig.  35b  reprecent  a 
segment  where  the  nonlinear  regression  failed  to  converge 
due  to  low  SNR  values.  The  circled  data  points  in  Pig.  33a 
are  measurements  with  low  wind  speeds  (2-3m/s)  and  small 
speed  differences  {c  O.Sm/s).  These  data  points  have  small 
speed  uncertainties  (  SVumii  0.2m/s)  which  cause  large 
values  of  6, •  The  boxed  data  points  in  Fig.  35s  exhibited 
very  broad  coherence  peaks  due  to  large  lateral  separation 
and  significant  changes  in  the  light  snow  inhomogeneities 
during  vertical  displacements.  Some  discrepancies  must  be 
expected  in  these  comparisons  since  lidar  measurements  are 
averaged  differently  than  conventional  wind  measurements. 

Messurement  of  wind  velocity  parameters  shown  in  this 
thesis  prove  thst  lidar  can  become  a  useful  tool  in  remote 
studies  of  the  boundary  layer.  The  reliability  in  lidar 
wind  measurements  can  be  increased  by  improvements  in  lidar 
instrumentation  to  decrease  the  noise  contributions  in  lidar 
returns  and  would  allow  wind  measxu'ementa  to  be  obtained 
\inder  s  greater  variety  of  atmospheric  conditions. 
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Appendix  A  Calculation  of  Spectra  and  Coherence 

A.l  Calculation  of  Spectra  and  Cross  Spectra  by  the  Fast 

Fourier  Transform 

The  Fast  Fourier  Transform  (KFT)  is  used  to  estimate 
the  power  and  cross  spectral  densities  of  discrete,  finite 
length  series  (Bingham  et  al.,  1967;  Cooley,  et  al.,  1970; 
and  Otneaa  and  Enochson,  197F}.  The  FFT  has  the  advantage 
over  previous  methods  of  spectral  estimation  in  computational 
speed,  especially  for  data  segment  lengths  that  are  integer 
powers  of  two. 

The  finite  site  of  the  data  segment  introduces  the 
problem  of  leakage  in  spectral  estimates.  Leakage  is  the 
effect  that  frequencies  or  wavenumbers  outside  a  particular 
interval  have  on  the  spectral  estimate  over  that  interval. 
Leakage  is  reduced  by  application  of  a  data  window  on  the 
segment.  Two  common  data  windows  are  the  Hann  window  and 
the  cosine  taper  (Otness  and  Enochson,  1972).  Both  windows 
decrease  the  side  lobe  characteristics  of  the  bandpsss.  The 
Cosine  taper  window  is  used  in  this  study  in  order  to 
balance  the  problems  of  statistical  stability  and  leakage  of 
the  spectral  estimates. 

The  cosine  taper  window  consists  of  tspering  ten  percent 
of  the  data  segment  length  on  each  end  with  a  cosine  bell 
given  by  the  formula  (Binghari,  et  al.,  1967), 
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D(x)* 


cos  (tfZx/L) 


~U2i  X  *  -/jL/10 
-UL/10  f  x<  UL/IO  A. 1.1 
JtL/10  <  X  t  L/2 


for  a  data  segment  extending  from  -L/2  to  L/2. 

The  data  window  is  applied  on  a  spatial  series^  V{x), 
in  the  following  manner. 


fUJ=f(x)  •  D(x) 


A. 1.2 


f(x)  is  the  modified  spatial  aeries. 

A 

The  raw  spectral  estimates  of  the  modified  series  f^, 
fg  •re  caicuiateo  irom  the  FFT  of  each  series  ^ f^  , 
respectiveljr. 

S^U)-  ^f^*  (k)  •  T  f^(k) 

S2(k)«  ?f2*  (k)  •  A. 1.3 

S  r  f^*{k)  •  '7' fjCk) 

Where  the  Fast  Fourier  Transform  is  defined  aa 

4«  I 

(  is  the  incremental  separation  of  rJ  data  points 
and  «  denotea  complex  conjugation.  $2  S^2 
respective  power  and  cross  spectral  estimates.  The  raw 
spectral  estimates  must  be  smoothed  in  order  to  decrease  the 
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P* 


mean  square  error  and  increase  the  statistical  signiricanee 
of  the  spectral  estimates.  Smoothinf;  can  be  accomplished  by 
a  i;eichted  average  over  adjacent  wavenarbers,  (Cooley,  et  al., 
1970),  or  by  time  averaging  over  a  number  of  modified  segments 
(Welsh,  19b7).  Time  averaging  was  chosen  for  this  study 
because  of  the  small  lengths  of  the  data  segments.  The 
smoothed  spectral  estimates  averaged  over  J  time  segments  are 


S^(k)  » 

j  f, 

S2(k)= 

A.l.u 

Si^k). 

•*  I 


Application  of  a  data  window  decreases  the  variance 
because  unequal  weighting  is  given  to  various  parts  of  the 
segment.  The  normalisation  factor  (d)  due  to  data  window 
application  is  (Otness  and  Enochson,  1972) 


Vs 

J  /  D^(x)  .  dJt 


A. 1.3 


Spectral  values  are  sealed  to  maintain  the  true  varianoa 


estimate. 

X/41 


.  J- 


A.l.b 


ASjSllt 
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A. 2  Calculation  of  the  Coherence 

The  coherence  between  two  time  or  spatial  series 
describe  their  relationship  at  various  frequencies  or  wave¬ 
numbers.  Calculation  of  the  coherence  between  two  series 
(subscripted  1,  2}  is  (Otness  and  tnochson,  1972} 


Sfp(k} 

Coh(k)=  - 


S^g(k} 


A. 2.1 


S^(k}  .  S^ik) 


where  S^(k),  S^lk)  are  the  smoothed  spectral  estimates  and 

smoothed  cross  spectral  estimate.  The  asterisk 
(e)  is  complex  conjugation.  The  calculation  of  smoothed 
spectral  estimates  are  described  in  Appendix  A.l.  Smoothing 
of  the  spectral  estimates  increases  the  degrees  of  freedom, 
therefore  decreasing  the  statistical  uncertaintj  in  spectral 
estimates  and  coherence  measurements.  The  value  of  the 
coherence  (P),  for  which  there  is  a  probability  (P)  of  a 
coherence  between  random  series  exceeding  this  value  is 
(Panofsky  and  Brier,  1966) 


r  .(l-p^/^***^*^}**  A. 2. 2 

where  edf  is  the  equivalent  number  of  degrees  of  freedom. 

For  a  chi-square  distribution,  this  number  is  twice  the 
number  of  averaging  segments.  In  this  study,  the  edf  is 
196  for  coherences  calculated  between  asinuth  angle 
(see  Fig.  2)  and  in  filter  estimations  (99  segments  sversgsd 


95 


over).  The  five  percent  probability  value  of  P  is 
.12.  Because  coherence  values  calculated  between  and 
(see  Fig.  2)  are  averaged  over  a  sr.all  nunber  of  segments 
(SO  )t  the  five  percent  probability  value  of  P  is  .17. 
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Appendix  B  Calculation  of  Hii^nal  Spectrun  to  foise  Cpcctru>n 
Ratio,  Sicnal  to  .''oise  Ratio  (SNR)  and  Optimal 
Linear  Filter. 

An  estimation  of  tne  ratio  of  the  signal  spectrun  to  the 
noise  spectrum  can  be  obtained  from  the  coherence  between 
different  spatial  series  separated  in  time.  The  spectra  of 
the  series  are  functions  of  wavenumber  and  time  lag  only.  The 
smoothed  spectral  estimate  composed  of  signal  and  noise  is 
(Weinstein  and  Zubokov,  1962) 

Tlk,  a  t)»  S^(k,  At)  ♦'s^{k,  At)  ♦  4t)  B1 

where  S(k,  At)  is  the  smoothed  spectral  estimate  of  the 
spatial  series,  S^(k,  At)  is  the  signal  spectrun,  S^(k,  At) 
is  the  noise  spectrun,  and  S  (k,  at)  is  the  "interference" 

8il 

spectmaa  of  the  signal  and  noise.  The  interference  tern  is 
zero  for  no  correlation  between  the  signal  and  noise.  Normal, 
random  noise  is  not  correlated  between  spatial  series  at 
different  times.  The  smoothed  cross  spectrxm  between  differ¬ 
ent,  spatial  aeries  is  the  cross  spectrum  of  the  signal  for 
normal,  random  noise,  i.e., 

a  t)- 

The  coherence  between  the  two  series  as  defined  in 
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Appendix  A is 

r  1  IL  n,i  ^a*?^**'*^^ 

Cohik,  dt)  *  :c - r: - 

S^(K,0)  •  S^lk,  it  ) 


b3 


where  «  denotes  complex  conjugation  and  the  subscripts  1,  2 
refer  to  the  two  separate  series.  In  the  limit  as  the  time 
lag  (^t)  approaches  zero*  the  cross  spectral  estimates 
limiting  value  is  tne  signal  spectrum.  The  individual 
spectra  of  each  series  approach  the  smoothed  spectral  estimate 
for  zero  time  lag, 

Lim  S.,p(k,  A  t)  =  S,,U,0) 

^  ^  B4 

Lim  Sp(k,  At)  ■  S^U.C) 

a  ^0 


The  limit  of  the  coherence  as  At  approaches  zero  is  there* 
fore  I 


Lim  Coh(k,  At)  « 
A  t*0 


S^l(k,0)  .  'S^jtk,©) 
S^U.O)  .  Sjik.O) 

S,^(k,0) 


2 


Se^ik.O)  ♦  S„iU.O) 


B5 


Solving  for  the  ratio  of  the  signal  spectrum  to  the 


noise  spectrum  yields 


Sgi(k,0) 


r Lim  Coh(k,  At)( 

tt-*0  J 


■  *s 


-1 


-1 
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Filtering  •  series  suppres&es  the  noise  component  in 
order  to  reproduce  the  signal  with  the  least  amount  of 
error.  The  filtered  series  ia  the  convolution  of  a  linear 

filter  (h)  with  the  original  series  (f)t 

•» 

Hx)  «  J^h(x^x*)  •  f(x*)  •  dx' 


The  optimal  linear  filter  that  minimises  the  quantity 


■  ^f(x)-f(x)J  ^  (overbar  denotes  ensemble  averages) 

is  calculated  from  the  knowledge  of  the  ratio  of  the  signal 
spectrum  to  the  noise  spectrum.  The  transfer  function  (H(k)) 
of  the  optimal  linear  filter  is  shown  by  Wainstein  and 
Zubakov  (1962)  to  be 

S^l(k,0) 

U(k)  »  -^3 - - -  B7 

s'^l(k.O)  ♦ 

where  ll(k)  is  defined  as  the  Fourier  transform  of  the 
optimal  linear  filter 

H(k)  «  / h(x*)  •  exp(-i  2  W  k  x»)  •  dx'»  Bti 

The  application  of  a  filter  changes  the  spectral 
eetimates  of  the  aeries.  The  Fourier  transform  of  the 
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filtered  series  is  the  product  of  the  transforms 

of  the  orlp.inol  series  and  the  filter  (K(k)j. 

0-'  ^ 
i  •  H(k) 

The  filtered  spectrum  is  the  conjugate  multipli¬ 

cation  of  the  Fourier  transform  with  itself. 


-  (H(k))^-  S(k) 


B9 


The  variance  of  the  series  is  the  integral  over  ell 
wavenumbers  of  the  spectrum. 

■  f  S[k)  •  dk 
o 

The  variance  of  the  filtered  series  is  thereforet 

r- 

(T/  «  /sik)-  dk  =  /(H(k)f  .  S(k)  •  dk 

•  ®  • 

For  white,  gaussian  noise,  the  noise  spectrum  is 
constant  over  the  bandpass  of  the  filter.  The  noise 
variance  of  the  filtered  aeries  is  given  byt 

-  /s' (k)  .  dk  «  T  /(HCk))^.  dk 

an  n  n  / 

S  being  the  constant  value  of  the  noise  spectrum, 
n 

The  weighted  integration  of  the  signal  spectrtBs  to 
noise  spectrum  ratio  by  the  square  of  the  transfer  function 


iTii 


is  the  retio  of  the  filtered  sicnal 

variance  (SNR). 

variance 

to  the  noise 

S^{k) 

SNR  «  J  (K{k)f  •  -  .  die 

*  “n<“> 

o 

.  /(H(k))^  .  [Htk)-*  -ij 

dk 

BIO 

'mr 
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Appendix  C  Derivation  of  an  Analytical  Form  for  Coherence 
Decay 

This  analysis  shows  that  the  mean  and  rms  wind  speed 
can  be  recovered  from  the  decay  of  coherence  between  lidar 
measurements  of  aerosol  inhoroogeneities  separated  in  time 
and  space.  These  measurements  are  assumed  to  extend  to 
infinity  in  the  presence  of  no  attenuation  in  order  to 
describe  a  simple  analytical  model  for  this  decay.  For  the 
purpose  of  this  analysis,  the  velocity  distribution  is 
assumed  to  be  homogeneous  and  isotropic  in  both  time  and 
space.  Tbe  aerosol  inhomogeneities  are  assumed  to  be 
described  as  three  dimensional  Gaussian  structures  being 
adveoted  by  the  wind. 

The  functional  form  of  the  inhomogeneities  is 

£>f[- k/i-  01 

where  t  refers  to  time  and  the  positions  are  described  in 
terms  of  generalised  coordinates,  *s.  The  initial 
position  of  the  centroid  of  the  inhomogeneity  in  the  ith 
direction  is  q^Q  with  a  rms  width  of  the  structure, C]- . 

The  notion  of  the  inhonogeneity  az*e  the  generalised 
velocities  in  each  direction,  q^. 

The  lag-cross  covariance  function,  C^2> 


between  two 
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liaor  |<rofilc3  ueparutcd  in  cpace  and  time  is  used  to  infer 
notion.  The  j/eomctry  of  the  experiment  is  shown  in  Kip,.  Cl 
for  covariance  In  the  direction.  One  lidar  profile  has 
q^  coordinates,  The  second  profile  is  separated 

from  the  first  in  tine  by  and  in  the  q^,  q^  directions 
by  The  lag-cross  covariance  averaged  over  all 

possible  initial  positions  and  time  is 

'■  .(  'a-  .1  %  L  ■(. 


The  variances  (C^,  C^)  of  the  profile  along  each  line  are 
averaged  in  a  similar  manner.  _ 

c,  ■  s  •<*  f  3 

o  •  as  '  •  A»  •  •'v 

C3. 

c.  s"*  /■%.  r^i. 

The  integrals  above  can  be  reduced  by  substitution  to  the 


form 


Ux  - 


sines  the  integrals  are  Gaussian  and  products  of  Gaussiana 


Plgur*  C«li  0«o««try  of  erots  corrolotion  in  qg  for  two  spatinl 
•orioo  bj  qj,  q3.  Qj,  !•  tho  log  v.luo  in  q^. 
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The  lar,-cross  correlation  coefficient,  ^  ,  defined  as 


Ci« 


is  after  all  integrations. 


c;>b 


The  ensemble  averaged  lag»cross  correlation  coefficient, 
over  a  normal,  three-dimensional  velocity  distribution. 


C(q)  is 


C6 


a(qj  is  assumed  to  have  the  form 

TTf^wO 


6ft  ^ 


C7 


where  ia  the  mean  of  the  velocity  in  the  ith  direction 
and is  the  rms  wind  speed.  The  integrands  are  iuuependeitt 
in  each  q^.  The  ensemble  averaged  lag-cross  correlation 
coefficient  is  the  product  of  the  separate  integrals 


u.iti 


Cda 
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The  integral  is  reduced  by  substitution  to  the  form  of 


F.q.  CJu  resulting  in 


Ctlb 


This  equation  Is  similar  to  one  used  by  Kunkel  (1978), 

To  simplify  the  functional  form,  the  aerosol  inhomopene- 
itics  arc  assumed  to  he  isotropic  (CT-  =  ^  =  )  and 

measurements  are  made  in  the  q^  plane  (4  q^=0).  The 
final  form  of  the  ensemble  averaged,  lag-cross  correlation 
coefficient  used  in  this  analysis  is  ^ 

Where  »  Ax,  Aq^  “Ay,  “  u,  U2  =  v,  :  O 

The  Fourier  transform  of  the  lag-cross  correlation 
coefficient  indy  is  the  lag-cross  spectrum,  ,ax,u,v,at) . 

The  integration  over  spatial  lag,  Ay,  is 

/  <rvr 

Which  can  be  reduced  to  the  sane  form  as  Eq.  Ck,  The 
resulting  lag-cross  spectrum  is 
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*4 


CIO 


The  coherence  function,  Coh,  ia  proportional  to  the 
aquare  of  the  absolute  magnitude  of  the  lag-cross  spectrum; 

(jc,»,c,'4«‘y 


The  value  of  the  coherence  for  the  quantity  #ix-uat>0) 
and  C”8«0,  ia  identically  one.  Therefore,  the  ratio 

^  ^J.  yjV<5*4<*J 


Cll 


ia  the  decay  of  the  coherence  with  spatial  and  temporal 
aeparationa  due  to  velocity  fluctuations.  Note  that  for 
■€,  ths  coherence  is  Gaussian,  as  expected  from  the 
assumptions  of  the  model. 
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Appendix  D  hrror  Arialysis  of  Pilot  Balloon  and  Tower 
Anenoineter  n’ind  ..easurencnts 

D.l  errors  in  Pilot  Balloon  lit uaurir.cnta  lAie  to  Position 
Errors 


Pilot  balloon  winds  are  derived  from  position  changes 
of  the  balloon  in  small  tine  intervals*  therefore  the  errors 
in  determining  the  position  of  the  balloon  are  translated 
into  speed  errors.  Tne  geometry  of  the  balloon's  position 
in  spherical  coordinates  is  shown  in  Fig.  D>1.  For  radar 
tracking  of  the  balloon*  the  slant  range*  H*  azimuth  and 
elevation  angles*/  fi  ttr«  measured.  The  horizontal  compon¬ 


ents  of  the  pibal  are*  (Middleton  and  Spilhaus*  1953) 


xsHcos(  ^)8ini/  ) 

yahcosi  O  )C08(/  ) 


01 


The  differential  element  of  each  component  is 

dxseosl^  l8in(#  )dP-K(sin(  9;8in(^  )dd  -cos(^  )cos(^  )d^) 
dy«cos(&  )eos(/  )dH-H(8in(  #  )cos(0  )d*  ♦coa(9)8in(^  )d/) 

02 

The  distance  tha  balloon  has  traveled  in  a  time  St  for  the 
X  component  is 

(u«$u)-^t»(x-»ix)/^  -(x^Sx)l^ 

replacing  exact  differentials  by  the  small  variations* 

The  noximun  velocity  error  in  the  x  component  is 


109 

Su... (/ 

-  vcM»d.s*^t<(  /  4. 

Jd 

—  /  C*/»  ©.  o-  f  ^  C**©,  U*0,  / 

* 

The  trigonometric  values  are  in  absolute  value  less  than 
one.  Therefore,  the  majtiinura  error  in  the  x  component  of 
velocity  is 

s<  %t  s<  ' 

where 

A  simila**  expression  can  be  derived  for  the  y  component 
of  the  velocity.  The  maximum  error  for  the  y  component  is 


fc’atimation  of  the  errors  in  terms  of  speed  and  aximuth  can 
be  obtained  in  the  same  manner  as  lidar  error  estimations. 
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D.2  Statistical  oncertainty  of  Tower  Ancnometer  Keasurenents 
Uncertainties  in  neasurenents  of  mean  quantities  by  an 
antmo.T.eter  depends  on  the  nunUer  of  independent  measurements 
made.  Keasurenents  of  this  Kind  are  time  averaf^ed  at  a 
single  i>olnt.  'itie  Eulerian  time  scale,  2^,  measures  the 
time  period  of  statistically  independent  samples  to  drift 
past  the  anemometer.  For  stationary  flows,  the  Eulerian 
time  scale  is  defined  as 


f  D.2.1 

o 

where  f  is  the  auto  correlation  function  and  dt  is  the  time 
lag.  The  number  of  independent  samples  in  a  time  period  of 

length  *7” la 


D.2. 2 


The  uncertainty,  C~  ,  in  determining  the  mean  value  of  a 
parent  population  with  a  standard  deviation, C*  and 
independent  neasurenents  is  (Bevington,  1969) 


G“  =  Cr/^  d.2, 3 

The  uncertainty  in  neasurenents  of  mean  quantities  by  an 


anemooieter  ia 


O'  *  ^ 
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